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ABSTRACT 

We are carrying out sensitive X-ray observations with Chandra and XMM of type II quasars selected 
from the Sloan Digital Sky Survey based on their optical emission line properties. In this paper 
we present the first results of our program. We present observations of four objects at redshifts 
0.4 < z < 0.8 and an analysis of the archival data for four additional objects in the same redshift 
range. Six of the eight were detected in X-rays; five of them have sufficient signal to derive spectral 
information. All of the detected sources have intrinsic luminosities L2-10 keV > 5 x 10''^ erg s~^; 
the five with sufficient counts for spectral fitting show evidence for significant absorption {Nh ^ a 
fewxlO^^ cm^^). At least three of the objects likely have Nh > 10^'^ cm^^; some may be Compton- 
thick {Nh > 10^^ cm~^). In the five objects for which we could fit spectra, the slopes tend to 
be significantly flatter than is typically observed in AGN; it is possible that this is due either to 
reprocessing of the nuclear emission or to a line of sight that passes through patchy absorption. 
Subject headings: galaxies: active — quasars: general — X-rays 



1. INTRODUCTION 

Type II Active Galactic Nuclei (AGN) are defined as 
those AGN lacking broad optical emission lines. In many 
cases, particularly at low luminosity (i.e., Seyfert galax- 
ies), the absence of broad optical emission lines has been 
shown to be due to obscuration, i.e., the central regions 
are blocked from the observer by a high column den- 
sity of gas and dust. In the standard unification model, 
the absorbing material is toroidal, resulting in highly 
anisotropic optical, UV and X-ray emission, while the 
strong narrow emission lines are produced above and be- 
low the obscuring torus in matter illuminated by the cen- 
tral engine. The detection of broad emission lines in po- 
larized, and hence scattered, light shows that at least 
in some cases the lack of broad lines is due to orienta- 
tion effects ( Antonucci .1993.') . Although low-luminosity 
obscured AGN (type II Seyfert galaxies) are abundant 
in th e local universe l|Kauffmann et al.ll 2004l : iHao et all 
120051^ . until recen tlv onlv a few ty pe II quasar candidates 
were known (e.g., Klcin mann et al.lll988i) . 

A key question is therefore to what extent the stan- 
dard model applies to type II quasars. This is the fourth 
in a series of papers using a sample selected on the basis 
of Sloan Digital Sky Survey (SDSS) data to study the 
numbers and properties of high-luminosity type II AGN, 
type II quasars. Our method is based on the supposition 
(as posited in the standard model) that most of the nar- 
row optical line emission originates outside the obscuring 
torus, and is therefore expected to be roughly isotropic. 
We have selected several hundred type II quasar candi- 
dates from the spectroscopic database of the S DSS based 
on th eir narrow emission line properties ifZakamska et all 
I200.'U hereafter Paper I). We are now conducting sen- 
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sitive multi-wavelength observations of these objects to 
determine the range of their properties and to under- 
stand how they compare to type II AGN in general. In 
this paper we present the first res ults of our Chandra 
ijWeisskopf et all 1200^ and XMM l|Jansen et alJ l200l|l 
programs. Our main goal is the determination of the 
X-ray spectral properties of this sample: their X-ray lu- 
minosities, continuum spectral slopes, and line-of-sight 
obscuration. 

Other methods have been used recently to search for 
type II AGN. Deep hard X-ray (i.e., E > 2 keV) surveys 
efficiently detect AGN candidates, including potentially 
obscured ones. For example, the X-ray background in 
the 1-6 keV band has bee n largely resolved into AGN 
at redshifts around z ~ 1 ijHornschemeier et al.ll200ll) . 
However, while many of these X-ray sources show only 
narrow emission lines in their optical spectra, a signifi- 
cant number of X-ray sources show broad emission lines 
or no emission lines whatsoever (e.g., Bargcr et al. 200^ 
Szokolv et alJ I2p04t^arger et all miOh: .Heckman et all 
20051 and lMatt)l20CI^ for an earlier review). In addition, 
many X-ray selected AGN candidates are very faint in 
the optical (i.e., R > 24). Consequently, it is difficult 
even to obtain complete redshift information for samples 
selected in this manner, much less a complete classifica- 
tion of their optical spectroscopic types. 

Alternatively Type-2 AGN are often detected in 
large-area hard X-ray surveys iD ella Ceca et al. 2002 
iCaccianiga et al.l 120041: iRevnivtsev et alJ 12004) and X- 
ray follow-u ps to very hard (E > 10 keV) X-ray all- 
sky surveys ijSazonov et al.ll2003) . These methods are 
promising since the detected objects tend to be (or are 
selected to be) bright in X-rays and often also in the opti- 
cal, allowing for spectroscopic classification and redshift 
determination. However, while several Type-2 quasars 
have been found (i.e., with Lx > lO'*'* ergs s~^), these 
sources tend to be mostly Seyfert 2 galaxies at low red- 
shifts (particularly in the very hard X-ray surveys since 
the high background of non-imaging detectors results in 
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high flux hmits). iSevergnini et al.1 l)2005j) discuss a re- 
lated approach in which a field with known Extremely 
Red Objects (EROs) was surveyed with XMM-Newton, 
and several EROs were found to harbor obscurred AGN 
and quasars. Likewise, a considerable effort is be ing ap- 
plied to use mid-IR colors to select AG N ((Haas et alJ 
I200I iLacv eral]l200l iStern et al.ll200^ . but the com- 
pleteness and efficiency of the selection of type II AGN 
by this m ethod has not yet been established. Along 
these lines. iMartinez-Sansigre et al.l l|2005|) selected type 
II quasars at z ~ 2 using the Spitzer First Look Survey 
data along with radio fluxes to find type II AGN and 
exclude starburst galaxies. This study suggests that the 
type II AGN fraction is ~ 50% at z ~ 2, but the un- 
certainties are large due to small number statistics and 
model-dependent assumptions. 

In Section [21 we review in detail the properties of the 
SDSS type II AGN sample. Section |3| describes our ob- 
servations and data reduction, Section 0] describes the 
X-ray spectral analysis. Section El gives details on indi- 
vidual objects. We discuss our results in Sectional fol- 
lowed by a brief summary in Section 7. An h — 0.7, 
flrn = 0.3, r^A = 0.7 cosmology is assumed through- 
out. We frequently refer to the [OIII]A5007A opti- 
cal emission line as simply [OIII]. Objects are identi- 
fied by their J2000 coordinates in Table El (e.g., SDSS 
J084234. 94+362503.1) and shortened to hhmm+ddmm 
notation elsewhere (SDSS J0842-h3625). 

2. SAMPLE DESCRIPTION 

The SDSS llYork et all l2f)f)0t ^toTighton et a,].l I200I 
lAbazaiian et alJl200 3. 2004. 20051) is an ongoing optical 
survey to image about 10,000 deg^ and obtain spectra 
of about 10^ galaxies and 10^ quasars. Using the spec- 
troscopic database of the SDSS as of July 2002 (about 
4 X 10'°^ spectra), we identified 291 objects in the redshift 
range 0.3 < z < 0.8 having strong narrow emission lines 
with high-ionization line ratios characteristic of an un- 
derlying AGN continuum (Paper I). Neither broad emis- 
sion lines nor the strong ionizing UV continuum are seen 
in these spectra, suggesting that the central engine is 
blocked from the observer. All these objects are there- 
fore type II AGN candidates based on their optical prop- 
erties. 

Since the narrow lines are presumably illuminated by 
the central engine similarly in type I and type II AGN, 
narrow line luminosities can serve as a proxy for the nu- 
clear luminosity. We used the [OIIIJA5007A emission 
line, which is present in all of the SDSS spectra of our 
objects, and one of the strongest optical emission lines 
in type II AGN l|Oster brock 1989). About 50% of the 
objects in our sample (130/291) have [OIII]A5007A lu- 
minosities in excess of 3 x 10^ Lq, similar to those of lu- 
minous unobscured quasars (—23 > Mb > —27). Based 
on this somewhat arbitrary criterion we classified objects 
with log(i[OIII]/Lo) > 8.5 as type II quasars; those 
with lower luminosities we designated type II Seyfert 
galaxies (Paper I). 

Using data from the IRAS all-sky survey we found 
that the mid-IR luminosities of the objects in our sam- 
ple reach a fewxlO**^ erg s^^ (Zakaniskact al. 2004), 
and therefore we directly confirmed the luminous na- 
ture of these objects. We have conducted spcctropo- 
larinietry of a subsample of the most luminous objects 



ijZakamska et al.l|2005j) and have shown that they harbor 
luminous blue broad-line quasars in their centers. These 
results are consistent with the basic orientation-based 
unification model of toroida l obscuration and off-plane 
scattering ( Antonucci 1993), implying that the model 
can be extended to include at least some high-luminosity 
AGN. 

We used the ROSAT All-Sky Survey (RASS; see Vo- 
ges et al. 1999 and references therein) data to investi- 

fte the X-ray properties of the objects in our sample 
akamska et al... 20041 . We found that type II AGN are 
about 10 times less likely to be soft X-ray sources than 
are type I AGN with the same redshift and [OIII] lu- 
minosity. Furthermore, an unexpectedly high fraction 
(50%) of the few objects with counterparts in RASS 
are radio-loud, in which cases the soft X-ray emission 
may be associa ted with jet activity rather than with the 
central engine ljUrrvl 1200^ . These findings indicate ei- 
ther that type II AGN are significantly undcrluminous 
in X-rays compared to type I AGN, or, more likely, that 
the soft X-ray emission is absorbed by intervening mate- 
rial in type II AGN, consistent with the standard uni- 
fication model. Similar conclusions were obtained by 
Vignafi et al. (2004), who investigated pointed ROSAT 
observations of the fields of 16 objects from our sample 
and found that the majority of objects are undetected 
in soft X-rays even in exposures much deeper than those 
of RASS. Our program of sensitive X-ray observations 
reaching well into the 2—10 keV band is intended to pro- 
vide better understanding of the low ROSAT detection 
rate of type II AGN, and to measure their absorbing col- 
umn densities and the X-ray luminosities independently. 

3. OBSERVATIONS AND DATA REDUCTION 

Of the eight objects presented in this paper, two were 
observed as part of our C/iondra proposal (ID 05701043), 
and two were observed as part of our XMM proposal 
(ID 020434). In addition, we checked all archival point- 
ings of Chandra and XMM within 15' of SDSS type II 
AGN^ and found that an additional four objects from 
Paper I were observed serendipitously. The list of the 
objects, observation IDs, exposure times and dates of 
the observations are given in Table 12 Of the four ob- 
jects observed serendipitously, two (SDSS J0115+0015 
and SDSS J0243+0006) do not meet our quasar luminos- 
ity criterion (L[OIII]> 3 x IO^Lq), so we classify them 
as Seyfert II galaxies. All eight objects in this paper are 
radio-quiet, as determined by the ratio of their radio and 
[OIII] luminosities (Zakamska et al. 200 4). 

We used XAssist l|Ptak fc Gr iffiths 20Q^ for initial 
processing of all data. XAssist works similarly for both 
Chandra and XMM-Newton data, with most operations 
performed by the software packages CIAO (version 3.2.1 
and CALDB 3.0.3) and XMMSAS (version 6.1.0). For 
each field, the data were reprocessed to take advantage 
of the latest calibrations, sources were detected, and the 
exposure was trimmed to remove background flares. The 
size of each source on the detector was estimated in or- 
der to determine appropriate source extraction regions, 
typically ~ 2" regions (Chandra) or ^ 16 — 18" regions 
(XMM-Newton) for on-axis point sources. For off-axis 

^ using the High-Energy Astrophysics Science Archive 
(HEASARC), http:/ /heasarc. gsfc.nasa.gov 



3 



observations we used larger source extraction regions, as 
described below for each individual case. Spectral re- 
sponses were calculated for the source spectra, and back- 
ground spectra were extracted from annuli centered on 
the sources, with interloping sources excluded from the 
background regions. 

4. RESULTS OF SPECTRAL ANALYSIS 

Of the eight objects in our sample, six were detected in 
X-rays, and all are consistent with being point sources. 
We defined "detection" as having enough counts in excess 
of the expected background that the Poisson probability 
of a spurious source at that location was < 1 x 10~^ 
(corresponding to a ^ 5ct criterion in Gaussian statis- 
tics). Since we are investigating specific locations rather 
than sampling a large field, this limit is very conser- 
vative. However, as discussed below, in only one case 
(SDSS J0842-h3625) would a lower (yet reasonable) de- 
tection limit result in a detection, and in that case the 
source detection is complicated by the position of the 
source in the field. For the purpose of estimating er- 
ror bounds for the luminosity, we treat all three weak 
sources, whether detected or not, in the same fashion. 

Sufficient counts for spectral fitting were detected from 
five objects (with more than 300 counts; see Table 2). 
We modeled the spectra with a power-law continuum ab- 
sorbed by the Galactic column density and an absorber 
at the redshift of the source, resulting in three free pa- 
rameters (the column density Nh, X-ray photon index F: 
dN/dE cx E-'^, and the power-law normalization). The 
relation between HI column density and opacity assumed 
solar abund ances. Galactic absorption was derived from 
the HI map l|Dickev fc Lockmanlll990j) using the nh tool 
provided by the HE AS ARC. The Galactic neutral hydro- 
gen column density does not exceed 4.8 x 10^" cm"^ for 
any of the objects in our sample. 

A simple absorbed power-law model is generally inac- 
curate for very high column densities of the absorber 
i^H ^ 10^'^ cm~^) because electron scattering sig- 
nificantly modifies the observed spectr um. We there- 
fore used the approximate model by lYaqooIj |1997f) 
that allows for scattering ("plcabs" in xspec). When 
the Compton optical depth is significant, the flux that 
emerges in any particular direct ion is strongly dependent 
on the scattercr's geometry fKrol ik. Madau. fc Zvckil 
[W94; GhiscUini, Haardt & Matt 1994); the plcabs model 
assumes spherical symmetry, which is not necessarily ap- 
propriate (we expect that the obscuration is toroidal). 
Our hope is that this model provides a slightly better es- 
timate, particularly for Compton depths that are order 
unity or less, than ignoring scattering altogether. For 
low optical depths, this model is equivalent to simple 
absorbed power-law models. 

The XMM-Newton spectra were binned in energy to 
20 counts per bin to allow the use of the statistic. 
For the two Chandra observations, the expected back- 
ground contribution to the source region is negligible 
(< 2 counts in both cases), and we used the C statis- 
tic l|Cashl IT979|) for spectral fits, without binning and 
subtracting the background. The parameter errors were 
computed at 90% significance for one interesting param- 
eter (Ax^ or C = 2.7). In the case of SDSS J0801-I-4412, 
which had only 40 detected photons, we fixed the power- 
law photon index at F = 1.7 (the typical value found 
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Fig. 1. — The Chandra ACIS spectrum of the serendipitously 
observed object SDSS J0115+0015, fit with an absorbed power-law 
with the absorber at the redshift of the source and with observed 
energies plotted. The residuals of the fit are shown in the bottom 
panel. 



in X-ray surveys for AGN fit with an absorbed power- 
law model; e.g., iNandra et all 120051 iTozzi et allboOSfl . 

Adopting values of F = 1.4 (the photon index of the 
X-ray background in the 0.5-10.0 keV bandpass), or a 
steeper slope such as F = 1.9 — 2.0, resulted in similar fit 
parameters. 

The spectral fits are shown in Figures d - 111 ^^nd Ta- 
ble 13 displays the parameters derived for the objects in 
our sample. The fits are statistically acceptable for all 
objects except SDSS J1641-f3858. For five of the eight, 
we were able to constrain the column density of absorp- 
tion, finding ~ 2 x 10^^ H cm~^ for four of them and 
roughly ten times that amount in the fifth case (albeit 
with poor statistics) . For four of the eight we were able 
to constrain the continuum slope, and all four are sig- 
nificantly harder than typical Type I AGN: three have 
F ~ 1.4, and the fourth has F ~ 0.5. After correction 
for light lost to absorption, the five spectral fits indicate 
2-10 keV luminosities between 0.8 and 7 x 10*** erg s~^. 

For two objects, SDSS J0243-I-0006 and 
SDSS J0842-I-3625, we regard our data as providing 
only upper limits on the X-ray flux because the Poisson 
probability of a spurious source at those locations was 
greater than 0.3%. Because SDSS J0243-h0006 was 
observed by XMM-Newton, each observation resulted in 
three separate upper limits . We therefore generalized 
the approach of iKraft et alJ 1)19911) . where upper fimits 
are derived for (single) observations in the Poisson 
statistical regime, to account for the different responses 
of the MOS and FN detectors as described in Appendix 
A. The results of this analysis are shown in Figure El 
which shows 99.7% confidence regions in the Lx — Nh 
plane for F assumed to be either 1.4 or 1.7. These 
regions are unbounded toward higher values of Nh since 
no fiux would be expected if Nh > 10^^ cm^^ (in the 
Compton-thick case, i.e., when the line-of-sigh t material 
is optically thick to Thompson scattering, iComastil 
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Fig. 2.— The XMM-Newton spectra for SDSS J0210-1001 fit 
with an absorbed power-law model. PN data are marked with 
circles, MOSl data are marked with triangles, and the MOS2 are 
marked with squares. The solid, dotted, and dashed lines in the 
upper-panel show the best-fitting model for the PN, MOSl and 
MOS2 data. 
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Fig. 3.— The Chandra spectra for SDSS J0801+4412 fit with 
an absorbed power-law model. The spectrum was binned here for 
display purposes only. 



1200^ . We plot the regions for both the absorbed and 
intrinsic X-ray luminosities. For comparison with the 
spectral fitting results, in Table |3| we also list the upper 
limits corresponding to fixing the column density at 
10^3 cm-2. 

We find it intriguing that the continuum slopes of the 
objects in our sample are consistently shallower than in 
type 1 AGN, e.g., F = 1.8—2.0 in citetrisa05, based on an 
XMM serendipitous survey of photometrically-selected 
SDSS quasars. It may be that this is because the spec- 
tral coni2lexity_often_se^^ of Seyfert II galax- 
ies l|Turner et alJll997t lRisalitill200a may be masked by 
our comparatively poor photon statistics. One possible 
complication is "reprocessing" , in which X-ray photons 
are reflected from optically-thick material located out of 
the line of sight (e.g., the accretion disk or the opposing 
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Fig. 4.— XMM-Newton spectra of SDSS J1226+0131 fit by an 
absorbed power-law model. Marks and lines are the same as in 
Figure!^ 
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Fig. 5.— XMM-Newton spectra of SDSS J1641+3858. Marks 
and lines are the same as in Figure 1^ 



side of the putative torus) or scattered by optically-thin, 
highly- ionized gas (e.g., the matter that creates the po- 
larization signature). Alternatively, the obscuring mate- 
rial may not uniformly cover the X-ray source. In both 
cases, the observed X-rays are absorbed by a range of col- 
umn densities, which tends to flatten the observed X-ray 
spectrum. 

To investigate this possibility, we have explored fit- 
ting a two-absorber model (both at the redshift of the 
source) to the two sources with the best photon statis- 
tics (SDSS J1226+0131 and SDSS J1641+3858). That 
is, in the limit of negligible Compton scattering {Nh <C 
10^"* cm~^), our model reduces to the form 

dN/dE = Aexp [-a{E)Nu,G.i] + z)]'^ x 
{/ exp [-~a{E [1 + z])NhM + (1 ^ /) exp [-a{E [1 -f z])Nh,2]} 

where E is the observed photon energy, a{E) is the ab- 
sorption cross section per H atom, / gives the fraction 
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Fig. 6. — The regions below and to the right of the curves 
are 99.7% confidence regions for X-ray luminosity (rest-frame 
2-10 keV) and column density for SDSS J0243+0006 (left) and 
SDSS J0842+3625 (right). X-ray luminosity uncorrected for ab- 
sorption is shown by dotted curves (F = 1.7) and dot-dash curves 
(F = 1.4). Absorption-corrected values are shown by solid curves 
(F = 1.7) and dashed curves (F = 1.4). For additional contrast, 
the F = 1.4 curves are colored red. Since these sources were not de- 
tected at high significance, these regions are generally unbounded 
from below (i.e., low luminosities are not excluded). In the case 
of SDSS J0842-I-3625, extremely low luminosity and high column 
density is formally excluded because the probability that this ob- 
ject was actually detected (99.5%) almost reached the confidence 
level used in this figure (99.7%). However, the range of column 
densities in question fall well into the regime where the approxi- 
mations of the model (plcabs) are questionable. We therefore treat 
these curves as showing only a luminosity upper limit. 



of the source obscured by column density Nh,i, while 
the remainder of the source is covered by matter with 
column density Nh.2- The two-absorber model is an ap- 
proximation to cither the case of partial covering or the 
case of some X-ray flux being scattered into the line of 
sight around a higher column-density absorber. The re- 
sults of the fits (i.e., the two column densities and the 
covering fraction / for the smaller column density) are 
shown in Table 01 For both objects, the best-fit column 
density on the line-of-sight of greater covering fraction in- 
creased, by a modest amount for SDSS J1641-I-3858, and 
by at least an order of magnitude for SDSS J1226-I-0131. 
In both cases the smaller column density covers < 10% 
of the source. As expected, the derived intrinsic con- 
tinuum slope grew significantly steeper in both cases 
(for SDSS J1641-|~3858, the resulting value of T is nom- 
inally larger than in the single-absorber case, but is not 
well-constrained). The results of fitting a two-absorber 
model should be treated with caution, because in the 
case of SDSS J1226-I-0131 the single-absorber model is 
statistically satisfactory (x^/dof= 28.1/25, which can be 
rejected at a confidence level of only 70%), so there is 
no requirement to introduce a more complex model in 
that case. However we note that these covering fractions 
(1-f, or > 90%) are simil ar to the values ob served in 
nearby Seyfert 2 galaxies (Turn er et al .11199^ . On the 
other hand, in the case of SDSS J1641-h3858, while the 
quality of fit is improved, it is still not entirely satisfac- 
tory (x^/dof= 99.8/75, which can be rejected at the 97% 
confidence level). We discuss this latter object in more 
detail in the following section. 

5. NOTES ON INDIVIDUAL OBJECTS 

SDSS J0115+0015: This AGN is located behind the 
cluster Abcll 168 and was covered serendipitously by a 
Chandra observation of this cluster, 1.75' off-axis. The 
X-ray source was detected with ~ 300 counts, with a neg- 
ligible background ('^ 2 counts in the source aperture of 



2"). Both its [OIII] luminosity (log(i[OIII]/io) = 8.14) 
and the intrinsic X-ray luminosity (6x 10'''^ erg s~^) place 
this object somewhat below the quasar luminosity cut-off 
(we adopted log(L[OIII]/Lo) > 8.5 as the optical defi- 
nition of a quasar from Pap er I, and Lx > 1 x 1 0'^'^ erg 
as the X-ray definition, iSzokolv et all 120041) . Nev- 
ertheless, this is still a luminous AGN and is therefore 
included in the sample. In addition to the fit to the un- 
binned spectrum, we also fit this spectrum after binning 
to 20 counts bin~^, which resulted in very similar fit pa- 
rameters. The single-absorber fit (Figure^ resulted in 
a column density of 2 x lO^'^ cm~^. 

SDSS J0210-1001: This object was observed as part 
of our XMM-Newton program. Its observed column den- 
sity is 2 X 10^^ cm^^ (the absorbed power-law fit is shown 
in Figure 121), however the observed spectral slope is very 
flat (F = 0.5). This may be a case in which reprocess- 
ing, patchy absorption or other effects are important (as 
described in Section but there are insufficient counts 
to justify any such model. 

SDSS J0243+0006: This object is within the field 
of view of two XMM-Newton observations of NGC 1068, 
located - 10' off-axis. Like SDSS J0115-f 0015, its [OIII] 
luminosity is not great enough for it to be labeled a 
quasar. It was not detected in any of the six exposures 
(using a source aperture of 21"), and the upper limits on 
the luminosity were obtained by combining the six expo- 
sures as discussed above and in the Appendix (see Figure 
EJ. Even summing the counts, the nominal Poisson prob- 
ability of no real source at this location was 15%. 

SDSS J0801+4412: This object was detected with 
~ 40 counts in a pointed Chandra observation from our 
program, too few to constrain the continuum slope. Less 
than 0.1 background counts are expected in the source 
aperture, so we used the C statistic with no binning and 
without background subtraction in our spectral fitting 
(Figure The resultant column density was '--^ 1 x 
10^'^ cm for assumed photon indices of either 1.4 or 
1.7. 

SDSS J0842+3625: This object has the sec- 
ond highest [OIII] luminosity, a fter IRAS 09104+4109, 
of all known type II quasars (|Kleinmam^_et_alJ [19^ 
Grawford & V anderriesti IT996I: iFranceschini et alJ 1200(1 
Iwasawa et alJ200lD a t redshifts z<l. Like IRAS 09104 
iHines fc Willsl 1199 3*). SDSS J0842-f3625 is 16% po- 
larized in the optical ( Zakamska et al. 2005). Unlike 
IRAS 09104, which is a double-lobed radio source, 
SDSS J0842-I-3625 has only a weak point-like radio coun- 
terpart and is radio-quiet, as determined by the radio-to- 
[OIII] ratio. The source was imaged by Chandra AGIS-I 
on two occasions, in Oct 1999 (observation ID 532, 8 ks 
exposure) and in Dec 2002 (observation ID 4217, 20 ks 
exposure) because it lies 6' away from the center of the 
cluster Abell 697. Unfortunately, in the latter exposure 
SDSS J0842+3625 happens to he between the AGIS-I de- 
tectors; we therefore discuss only the 8 ks observation. In 
that image, SDSS J0842+3625 is 6' off-axis, where the 
point-spread function of Chandra is rather larger than 
on-axis. We used a 4" extraction radius around the 
SDSS position, in which we found six counts. To esti- 
mate the background, we used two circular regions 25" 
in radius offset from the source position, from which we 
estimated there would be 1.0 ± 0.2 and 1.5 ± 0.2 back- 
ground counts in the source region. This variation, while 
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significant at less than the 2cr level, may be due to the 
cluster emission. The count-to-background ratio from 
the higher background estimate implies that the source is 
detected with 99.5% confidence (Kraft ct al. 1991). We 
therefore conservatively consider the observation to be 
an upper- limit ( see Figure IBIl . This source is described 
as undetected bv lVignali et al.l l)2004|) based on the same 
data. 

SDSS J1226-I-0131: This object was serendipitously 
covered by an XMM-Newton observation of an extra- 
galactic gas-rich system, HI1225-I-01. A significant frac- 
tion of the exposure (~ 50%) was affected by strong 
flaring, and the data collected during the flares were ex- 
cluded from the analysis. The single- absorber fit resulted 
in a column density of Nh = 2 x 10^^ cm~^ (Figure 

B . The spec t rum o f this object was also presented by 
ignali et al] l)2004j) . who obtained similar spectral pa- 
rameters. 

SDSS J1232+0206: This object was a target in our 
Chandra program. Six counts were detected within a 2" 
aperture at the position of the source and the background 
estimate in the source aperture is 0.3 counts. This means 
that the source was marginally detected (Poisson prob- 
ability of a false detection 1 x 10^^). Since we have 
only barely enough counts for a detection, we follow the 
procedure discussed in the Appendix. We find that the 
99.7% confidence contours (equivalent to 3a in Gaussian 
statistics) define a band in the luminosity-column density 
plane (Figured). Because it is detected, albeit weakly, 
we can exclude the very low luminosity and high col- 
umn density corner of this plane. We do not consider 
X-ray luminosities below 10**^ erg s~^ because starburst 
X-ray emission can dominate over low-luminosity AGN 
emission at this level (e.g., see Ptak et al. 1998). In ad- 
dition, we demonstrated in our earlier work (Zakamska 
et al. 2005) that scattering in the interstellar material in 
the host galaxy on kpc scales can be a significant effect, 
and around 1% of the intrinsic optical emission can thus 
reach the observer. If scattering is primarily by electrons, 
then the same scattering efficiency applies to the X-ray 
emission, and an apparent X-ray luminosity of 10"*^ erg 
can be dominated by the scattered light that has 
not been subject to circumnuclear obscuration. There- 
fore, for Nh S 10^^ cm^^, our upper limit effectively 
becomes Lx S 10^"^ org s^^. This object was listed by 
Zakamska et al. (2004) as having a potential counter- 
part in the RASS catalog, with a nominal offset between 
the SDSS position and the RASS position of 0.75'. Our 
Chandra data indicate that this is likely a mismatch, as 
there is a much brighter source in the field closer to the 
nominal RASS position. 

SDSS J1641+3858: This object was observed by 
XMM-Newton as part of our program, and ~ 1600 counts 
were detected. As mentioned above, a single absorber fit 
(Figure is statistically unacceptable. The addition of 
a second absorber resulted in a marginally acceptable 
fit. Residuals remained at an observed energy of ~ 4.0 
keV, which corresponds to 6.4 keV in the rest frame, the 
energy of the neutral Fe Ka emission line. We added 
a narrow (0.01 keV physical width) Gaussian to the fit, 
which reduced by ~ 8 for 2 additional parameters, 
but this does only slightly improves the confidence level 
of the fit. The rest-frame energy of the line was 6.69 
(6.60-6.80) keV and the equivalent width (EW) was 0.16 
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20 21 22 23 24 

log N„ (cm-2) 

Fig. 7. — 99.7% confidence regions for X-ray luminosity (rest- 
frame 2-10 keV) and column density for SDSS J1232+0206. X-ray 
luminosity uncorrected for absorption is shown by dotted curves 
(r = 1.7) and dot-dash curves (F = 1.4). Absorption-corrected 
values are shown by solid curves (F = 1.7) and dashed curves 
(F = 1.4). For additional contrast, the F = 1.4 curves are colored 
red. As discussed in the text, this source was detected at > 99.7% 
confidence and therefore the confidence region is bound from below 
(i.e., very low luminosities are excluded). 

(0.05-0.35) keV. This line energy is signicantly higher 
than 6.4 keV, which is expected from "neutral" Fe (less 
ionized than Fe XVII) and is cons isten t with He-like Fe- 
Ka emission. iNandra et al.l l)1997|) and lReeves fc Turned 
((2000) found evidence for a luminosity dependence in 
the Fe-K emission of AGN, possibly due to the puta- 
tive accretion di sk becoming ionized a t high accretion 
rate. However, I.Timenez-Bailon et al.l (|2005) found no 
correlation between Fe-K line energy and X-ray lumi- 
nosity in a sample of PG quasars. The observed EW 
is similar to values obse rved in Type-1 radio-quiet AGN 
(|R eeves fc Turnei<l2'000l) . suggesting that the continuum 
we are observing at 6 keV is not obscured. This fit is 
shown in Figure |H1 

6. DISCUSSION 

The primary goal of our study was to determine the X- 
ray spectral characteristics of a sample of type II quasars 
which were selected based on their narrow emission lines. 
Because our sources were all identified exclusively on 
the basis of optical properties, our sample should be 
free of any X-ray selection bias. The four in our pro- 
gram were chosen as those having the greatest [OIII] fiux; 
SDSS J0842-I-3625 would have been in our program but 
for the fact that it had already been observed in a clus- 
ter observation; the other three were observed serendip- 
itously. 

6.1. The column density distribution 

All five of our type II quasars with enough counts to 
constrain absorption showed column densities of at least 
^ 10^^ cm^^. As we will argue below on the basis of 
its X-ray/ [OIII] luminosity ratio, we believe it likely that 
SDSS J0842+3625 was undetectable because its obscura- 
tion is Compton-thick. Thus, while our sample is small, 
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Fig. 8.— The XMM-Newton spectra of SDSS J1641+3858 fit 
with a power-law plus two absorbers and a narrow Gaussian. The 
Gaussian line energy was consistent with a rest-frame energy of 
~ 6.7 keV. Note that this fit is marginally statistically acceptable 
(rejected at a confidence level of 93%) while the single-absorber 
power-law fit is clearly not (rejected at a confidence level of 99.5%). 



we have nonetheless shown it is very hkely that type II 
quasars are in general strongly absorbed. 

Our sample can be compared to X-ray observations of 
nearby Seyfert II galaxies, with the caveat that archival 
X-ray surveys of nearby AGN suffer from an X-ray se- 
lection bias because the AGN were usually observed at 
least in part on the basis of their X-ray flux. In addition, 
it should be kept in mind that at lower luminosities the 
soft X-ray band can be dominated by starburst emission. 
Several such surveys exist in the literature.^ 

In a n ASCA survey of Seyfert II galaxies. iTurner et al.l 
l)1997|) found somewhat steeper continua, but very simi- 
lar column densities, to those encountered in our sample: 
the mean power-law index and absorbing column in their 
study were T — 1.8 and Nh = 4 x 10^^ cm^^. However, 
often the absorbed power law fits were poor. Fitting their 
data with more complex models, such as partial covering 
with neutral or ionized material, a thermal plasma com- 
ponent (representing starburst emission below 2 keV), 
or refiection (off neutral, Gompton-thick material out of 
the line of sight), improved the fits. In the case of partial 
covering by neutral material, the mean photon index in- 
creased to 2.0 and the mean column density increased 
to 1 X 10^^ cm~^, although in a number of cases a dou- 
ble power-law model (in which the two power-laws were 
allowed to have different slopes and the absorption for 
one power-law was fixed at the Galactic value) provided 
b etter fits . 

iRisalitil l)2002() presented a BeppoSAX survey of 
Compton-thin Seyfert II galaxies, where the simplest 
model fit to the data was a double power-law model. 
The mean spectral index and column density (for the ab- 
sorbed power-law) from this analysis were F ~ 1.8 and 



Nh 10^3 cm"2^ similar to the ASCA resuhs. Note 

that the partial-covering fits to the ASGA data resulted 
in a slightly steeper slope (F ~ 2) than found here, prob- 
ably due to the inclusion of PDS data at higher energies, 
where reflection effects become significant. 

The previ ous work mos t near ly comparable to ours 
was that of iRisaliti et al.l l)1999|) . who studied the X- 
ray column density distribution of a sample of type II 
Seyfert galaxies, likewise selected on the basis of [QUI] 
flux. These authors found that 75% had column densities 
Nh > 10^'^ cm^^ and 50% were Gompton-thick. Thus, 
if anything, the evidence in hand suggests that while the 
incidence of absorption in type II Seyferts and quasars 
is similar, its thickness may be somewhat greater in the 
Seyferts. However the statistics are very limited and a 
larger sample is necessary to confirm the trend and to 
determine whether this contrast (if real) is primarily due 
to luminosity or redshift. 

6.2. Spectral slopes 

We were able to constrain the continuum slope in 
four objects and achieve an acceptable fit for three of 
the spectra. As we have already remarked, in general 
we find rather harder continuum slopes than are com- 
monly found in type I AGN, whether Seyfert galaxies or 
quasars. As mentioned in the previous subsection, our 
slopes are also harder than those found in previous stud- 
ies of type II Seyfert galaxies. 

This contrast in slopes may not be genuine. Studies of 
type II Seyfert galaxies with better photon statistics have 
shown that more complicated models than a power-law 
with a single abso rber are often requ ired to adequately 
describe the data l|Turner et alJll997(l . 

By fitting the two objects with the best statistics with 
more complex spectral models, we have seen indications 
that these hard spectra (and possibly also our mod- 
est column densities) are indeed artifacts of fitting low 
signal-to-noise data with a model simpler than is appro- 
priate. However, we do not have sufficient signal to prop- 
erly test this speculation. 

6.3. [OIII] emission as an indicator of X-ray luminosity 

Tables El and H list the Lx/L[OlU] ratios. The [OIII] 
luminosities of the objects in our sample are not cor- 
rected for possible reddening of the narrow-line region, 
since the Ha emission line is redshifted out of the op- 
tical range for all objects in our sample, and thus we 
are not able to compute the narrow-line Balmer decre- 
ment. These r atios may be compare d with those found 
in the study of'Heckma n et all l)2005j) , who analyzed the 
hard X-ray and [OIII] luminosities of nearby AGN. They 
found that when a sample of type II Seyferts is selected 
by [OIII] fiux, the mean log^oiLx / L[Om]) = 0.6 with 
standard deviation 1.1. All our detected objects (and the 
upper limit for one of our undetected objects) fall well 
within the extrapolation of this distribution to higher 
luminosities (Figure EJ. In this respect, our type II 
quasars fall nicely in line with the behavior of type II 
S eyfert galaxies. 

iHeckman et al.l l)2005j) found in addition that the mean 
X-ray/ [OIII] ratio of hard X-ray-selected type II Seyferts 
is rather larger than when the selection is based on [OIII] 
flux, and is then quite similar to the ratio for similarly- 
selected type I Seyferts: {log{L x / L[Olll])) = 2.2. The 
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Fig. 9. — 2-10 keV X-ray luminosity (in units of ergs s^^) plotted 
as a function of L[OIII] (in units of Lq). Seyfert I (empty dia- 
monds) and Seyfert II (filled diamonds) samples were taken from 
IHeckman eTaU l200l) . The type II quasars from this paper are 
shown with filled and empty circles showing the X-ray luminosity 
before and after correcting for absorption (as listed in Table 2), 
i.e., in both the Seyfert and our AGN samples, the empty symbols 
should represent intrinsic X-ray luminosities. The dashed and solid 
lines s how the mean log{L x / L[OlU]) values from Hcckman et al. 
12003) for the [OIII]-selected Seyfert I and II samples (1.6 and 
0.6 dex, respectively; neither Lx nor L[OIII] were corrected for 
absorption in these samples). The upper-limits are based on as- 
sumed column densities of 10^'^ cm~^, and the corresponding limits 
on intrinsic luminosities are a factor of ~ 2 higher than shown. 



most natural interpretation of this result is that many 
type II Seyferts identifiable through [OIII] emission are 
absent from hard X-ray-selected samples, presumably be- 
cause even hard X-rays are entirely blocked by Compton- 
thick obscuration. A similar conclusion was found by 
iBas sani ct al. (1999), where it was shown that the Fe-K 
emission line equivalent width in Seyfert 2 galaxies in- 
creases strongly with decreasing X-ray/[OIII] ratio, with 
Compton-thin AGN mostly exhibitng X-ray/ [OIII] > 1. 
This also suggests that the X-ray continuum is largely 
becoming absorbed when the observed X-ray/ [OIII] ra- 
tio is low. Because our upper bound on Lx/-Zj[OIII] 
for SDSS J0842-I-3625 (assuming a column density of 
10^'^ cm^^) is nearly two orders of magnitude below the 
typical ratio for hard X-ray selected AGN, we suggest 
that it may be such a case of Compton-thick obscura- 
tion. 

In fact, although our detected objects are roughly con- 
sistent with the distribution of Lx / L[Oni] ratios for 



[OIII]-selected type II AGN, they fall systematically be- 
low the mean relationship for X-ray selected AGN. This, 
too, suggests that we may be underestimating the true 
absorption even in these cases. 

7. SUMMARY 

Our results have in large part confirmed that the stan- 
dard unification model satisfactorily describes optically- 
selected type II quasars. Six of the eight are strong X- 
ray sources, and all five with enough counts to fit spectra 
show substantial absorbing columns. Outstanding ques- 
tions include: Is the contrast in continuum slopes that 
we find between typical type I quasars and these ob- 
jects real? Is the column density distribution of type II 
AGN a function of redshift or luminosity, as hinted by 
this small sample? Observations of larger samples with 
greater signal-to-noise ratios will be necessary to make 
progress. 
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APPENDIX 

MULTI-DETECTOR LUMINOSITY AND COLUMN DENSITY LIMITS 

As discussed in Kraft et al. (1991) and Van Dyk et al. (2001), Bayesian statistics are particularly useful for the 
interpretation of data with low numbers of counts. We begin by computing the Bayesian posterior probability for 
X-ray luminosity (Lx), column density (Nh), photon index P and background within the source aperture (B), based 
on flat (i.e., constant) priors on Lx and Nh- 

n 

P{Lx,NH,r, B,„r,\D) = CP(r)l[C{T,\S,{Lx,NH, P) + B,)C{Tb,,\B,A,/As) (Al) 
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where C{Ti\Si{Lx, Nh, T) + i?i..„) is the hkehhood of observing Ti total counts in the source aperture from observation 
i, given a count estimate Si + Bi, and £{TB,i\BiAb/As) is the likehhood of observing Tb^i counts in the background 
region (with area Ab) when the background estimate is B in the source region (with area As). The product is over the 
individual observations i, with each XMM detector being considered a separate observation since the responses differ. 
Both likelihoods arc Poisson distributions. Si is the number of counts expected given the spectral model parameters, 
which is computed for the response (calculated at the position of the source) using xspec. P{T) is the assumed prior 
distribution of photon indices, however here we fix F at the values of 1.7 and 1.4 (i.e., equivalent to a 5-function prior). 

We estimate Bi using an annulus around the source or a source- free region close to the source. C normalizes the pos- 
terior probability to unity. The background counts estimates -Bi..„ are nuisance parameters which can be marginalized 
from the posterior (hereafter dropping the \D notation for the posterior): 

P{Lx,NH,r)^ dBi dB2-- dBnP{Lx,NH,T,Bi„„)^T\ dB,P{Lx, Nh.T , B,) (A2) 
Jo Jo Jo j Jo 

The posterior can then be used to compute confidence contours in the Lx — Nh plane, which would be the smallest 
regions that encompass the desired confidence level. In practice, when we have an upper limit on the count rate, the 
most probable values will be near the highest allowed Nh and lowest allowed Lx, and the confidence contour will be 
unbounded. 
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TABLE 2 

SDSS TYPE II AGN OBSERVED WITH Chandra or XMM-Newton 





J2000 




Galactic N^j 






Observation 


Exposure 


Date 


off-axis 




coordinates 


(xlO^" cm-2) 


z 


Iog{L[OIII]/L0) 


ID 


(ks) 


mm/dd/yy 


angle (') 


SDSS 


J011522.194 


-001518.5 


3.4 


0.390 


8.14 


Chandra-mOi 


37.6 


11/01/02 


2.0 


SDSS 


J021047.01- 


■100152.9 


2.2 


0.540 


9.79 


XMM-0204340201 


9.7(P),11.6{M1),11.6(M2) 


01/12/04 




SDSS 


J024309.794 


-000640.3 


3.6 


0.414 


7.95 


JfMM-0111200101 


35.3{P),38.7(M1),35.6{M2) 


07/29/00 


10.8 














XMM-0111200201 


33.0(P),37.8(M1),34.9{M2) 


07/30/00 


10.4 


SDSS 


J080154.244 


-441234.0 


4.8 


0.556 


9.58 


Chandra-52A8 


9.8 


11/27/03 




SDSS 


J084234.944 


-362503.1 


3.4 


0.561 


10.10 


Chandra-532 


7.3 


10/21/99 


5.4 


SDSS 


J122656.484 


-013124.3 


1.8 


0.732 


9.66 


XMM-01 10990201 


9.0(P),15.3{M1),15.4(M2) 


06/23/01 


6.0 


SDSS 


J123215.814 


-020610.0 


1.8 


0.480 


9.69 


Chandra-A91l 


9.5 


04/20/05 




SDSS 


J164131.734 


-385840.9 


1.2 


0.596 


9.92 


XMM-0204340101 


12.2(P),16.8(M1),17.1{M2) 


08/20/04 





Note. — J2000 coord inates, rcdshifts and [OIII]A5007A luminosities are from Paper I. The Galactic absorption was derived from the HI map by 
iDickev & LockmanI fl990D using the nh tool provided by the HEASARC. For XMM-Newton observations, the exposure times are Ustcd separately for PN (P) 
and MOSl,2 {Ml, 2) instruments. The date of the observation is given in the next column. The off-axis angle is given in arcmin only for those four objects 
that were observed serendipitously. The remaining four objects were the primary targets of the observations. 



TABLE 3 

X-RAY SPECTRAL PROPERTIES OF SDSS TYPE II AGN 



ID 


Counts 


Nh 


r 




Lx 


Lx.i7itr 


Lx/i[OIII] 


Lx,i„tr/L[Om] 






(1022 cm-2) 






(10^4 erg s-i) 


(10^4 erg s^l) 




SDSS J0115+0015 


339 


2.0 (1.3-2.8) 


1.39 (1.09-1.73) 


471/523 


0.71 


0.82 


134 


155 


SDSS ,10210-1001 


300 


2.3 (0.4-5.5) 


0.46 (-0.07 - 0.72) 


9.1/9 


1.8 


2.0 


7.5 


8.3 


SDSS J0243+0006 




10 (fixed) 


1.7 (fixed) 


< 0.060 


< 0.10 


< 18 


< 30 


SDSS J0801+4412 


40 


16 (9.7-27) 


1.7 (fixed) 


221.9/524* 


0.87 


1.9 


6.0 


13. 


SDSS J0842+3625 


< 15 


10 (fixed) 


1.7 (fixed) 




< 0.47 


< 0.81 


< 1.0 


< 1.7 


SDSS J1226+0131 


574 


2.0 (1.4-2.7) 


1.41 (1.16-1.69) 


28.1/25 


4.0 


4.6 


23 


26 


SDSS J1232+0206 


5.7 


10 (fixed) 


1.7 (fixed) 


0.25 


0.45 


1.3 


2.4 


SDSS J1641+3858 


1624 


2.8 (2.2-3.4) 


1.33 (1.19-1.51) 


113.5/77 


5.7 


6.7 


18. 


21 



Note. — The counts given in the case of XMM data arc the net counts from the three detectors combined. In the ease of SDSS J0243+0006, no counts 
value is given since the upper limit on counts combined from disparate detectors is not meaningful. The errors on fit parameters Nh (hydrogen column 
density at the redshift of the source) and F (photon index, dN/dE oc E~^) are based on Ax^ = 2.7. Lx is the rest-frame 2—10 keV luminosity and Lx,intr 
is the rest-frame 2 — 10 keV luminosity after correction for absorption. 



Fitting was performed using the C statistic with unbinned data. 



TABLE 4 
Partial-Covering Fits 



ID 


Nh,i 


Nh,2 


r 


/ 


x2/dof 


Lx /Lx.intr 


Lx/L[01U]/Lx,intr/L[0UI] 




(1022 cm-2) 


(1022 cm-2) 






(10"" erg s-i) 


SDSS J1226+0131 


3.6 (3.1-5.1) 


87 (53-150) 


2.83 (1.96-3.81) 


0.091 


19.6/23 


4.0 / 40. 


23 / 230 


SDSS J1641+3858 


0.41 (0.10-1.8) 


4.6 (3.4-8.7) 


1.56 (1.39-1.78) 


0.071 


99.8/75 


5.7 / 7.6 


17.9 / 23.8 



Note. — Results of fitting a model consisting of two absorbers plus a power-law to the highest signal-to-noise ratio spectra. The quantity / 
is the fraction of fiux intercepted by the absorber with the lower column density. See text for details. 



